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Casimir effect mechanism of pairing between fermions in the vicinity of a magnetic 

quantum critical point 
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We consider two immobile spin 1/2 fermions in a two-dimensional magnetic system that is close 
to the 0(3) magnetic quantum critical point (QCP) which separates magnetically ordered and 
disordered phases. Focusing on the disordered phase in the vicinity of the QCP, we demonstrate 
that the criticality results in a strong long range attraction between the fermions, with potential 
V(r) oc —Ijr'^, V « 0.75, where r is separation between the fermions. The mechanism of the 
enhanced attraction is similar to Casimir effect and corresponds to multi-magnon exchange processes 
between the fermions. While we consider a model system, the problem is originally motivated by 
recent establishment of magnetic QCP in hole doped cuprates under the superconducting dome at 
doping of about 10%. We suggest the mechanism of magnetic critical enhancement of pairing in 
cuprates. 
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I. INTRODUCTION 

In the present paper we study interaction between 
fermions mediated by magnetic fluctuations in a vicin¬ 
ity of magnetic quantum critical point. To address this 
generic problem we consider a specific model of two holes 
injected into the bilayer antiferromagnet. The results 
presented below demonstrate that critical magnetic fluc¬ 
tuations lead to the long range Coulomb-like attraction 
between the holes. 

Our interest to this problem is motivated by cuprates. 
Lying at the center of the debate of high-Tc superconduc¬ 
tivity is whether it originates from a Fermi liquid or from 
a Mott insulator. Recent experimental data, including 
Angle-Resolved Photoemission Spectroscopy (ARPES) 
support Mott insulator scenario in underdoped cuprates 
and show transition from small to large Fermi surface 
in the hole doping range 0.1 < a; < 0.15, see Refsii^— 
Magnetic quantum oscillations (MQO) in underdoped 
YBa 2 Cu 306 -i-j/ also support the small pocket scenario^, 
in contrast to the large Fermi surface observed on the 
overdoped side^. Besides that, existence of hole pockets 
is consistent with the picture of dilute holes dressed by 
spin fluctuations, based on doping a Mott insulatori^ 

Optimally doped and overdoped cuprates do not have 
any static magnetic order. On the other hand, the un¬ 
derdoped cuprates possess a static incommensurate mag¬ 
netic order at zero temperature. A magnetic QCP sep¬ 
arating these two regions was predicted in Refi^ at dop¬ 
ing X ~ 0.1. In La 2 -xSra;Cu 04 the QCP is smeared 
out because of disorder. However, in YBa2Cu306-i-j/ the 
QCP is located experimentally with neutron scattering, 
nuclear magnetic resonance wipeout and muon spin rota¬ 
tion (/iS'i?) at doping x « 0.09 {y ~ 0.47)A"— At larger 
doping, after crossing the QCP the (quasi-) static mag¬ 
netic ordering vanishes and becomes fully dynamic. 

It is widely believed that superconducting pairing in 
cuprates is driven by a magnetic mechanism. The most 
common approach is based on the spin-fermion model 


in the frame of normal Fermi liquid picture (large Fermi 
surface). Within this approach electrons interact via ex¬ 
change of an antiferromagnetic (AF) fluctuation (param- 
agnon) The lightly doped AF Mott insulator approach, 
instead, necessarily implies small Fermi surface. In this 
case holes interact/pair via exchange of the Goldstone 
magnonJ^ Due to the strong on-site Hubbard repulsion 
both approaches result in the d-wave pairing of fermions. 

Magnetic criticality can significantly influence super¬ 
conducting pairing. This idea has been recently consid¬ 
ered by Wang and Chubuko\d^ in a context of electron 
doped cuprates. There are also some earlier works ref¬ 
erenced in Refji^. However, to the best of our knowl¬ 
edge all the previous works imply a normal liquid with 
large Fermi surface. This might be a reconstructed Fermi 
surface which emulates small hole pocketsi^, but still in 
essence this is a weak coupling normal Fermi liquid like 
approach. A large Fermi surface to a significant extent 
diminishes importance of the magnetic criticality for the 
pairing. 

In this work we consider two holes injected in the 2D 
“rigid” Mott insulator, so in essence our approach implies 
the small Fermi surface and therefore strong coupling 
limit. In this case influence of the magnetic criticality on 
the coupling between two fermions is the most dramatic. 
As the Mott insulator host we use the bilayer antiferro¬ 
magnet with magnetic fluctuations driven by the inter¬ 
layer coupling. We consider the bilayer model for the 
sake of performing a controlled calculation. However, 
we believe that conceptually our conclusions are equally 
applicable to the single-layer and multi-layer cuprates. 
The model presented here has only commensurate mag¬ 
netic ordering, so we put aside incommensurability in the 
cuprates. 

The model under consideration demonstrates spin- 
charge separation at the QCPj^^ It means delocalization 
of hole spin due to dressing by divergent magnon cloud. 
The effect of spin-charge separation points out to the 
nontriviality of the pairing problem. We are not familiar 
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with any models that consider pairing of two fermions, 
that incorporate physics of spin-charge separation. 

In order to probe the interaction between two fermions 
we consider spin fluctuations in the system, keeping the 
fermions to be immobile and spatially localized, just as 
magnetic impurities Because of the immobility, Fermi 
statistics of the impurities is not actually relevant for our 
results, however for shortness we will call the impurities 
as ’’fermions” hereafter. Our calculations show that the 
single magnon exchange is getting irrelevant close to the 
QCP. Instead, we obtain strong inter-fermion attraction 
in singlet and triplet spin channel due to Casimir effect.— 
Each of injected fermions (holes) builds up a ’’bag” of 
the quantum magnetic fluctuations. The fermions at¬ 
tract to each other, sharing common bag and reducing 
energy of the magnetic fluctuations inside of the bag. 
So-called ’’spin-bag” mechanism of attraction in the con¬ 
text of high-Tc superconductivity in cuprates is famil¬ 
iar from the works of J. Schrieffer et aZ.— , that however 
did not account for the magnetic criticality, neither small 
Fermi surface. The spin bag model has also conceptual 
similarity to QCD bag models for nucleon binding such 
as MIT— and chiral bags^ that have being extensively 
studied from I970’s up to now. 

The structure of the paper is following. In the Section 
m we introduce bilayer J — J± antiferromagnet, which 
is simple but instructive model and contains all essen¬ 
tial physics of magnetic criticality. In the Section III Al 
we characterize magnetic quantum critical point driven 
by interlayer coupling J_l/J and describe magnon exci¬ 
tations for undoped AF in disordered phase in the frame 
of spin-bond mean held theory. Next, in the Section Hi Bl 
we move to the hole-doped J — J± model and show how 
holes interact with magnons. In the Section ITTTl which is 
the main content of the paper, we consider hole-hole pair¬ 
ing problem at the QCP and show that pairing can not 
be described in terms of one-magnon exchange. In Sec- 
tion llll Al we develop effective theory for Casimir interac¬ 
tion of the fermions, considering double-fermion ’’atom” 
which can be either in singlet or triplet state. In Section 
IIII Bl we present results of solution to Dyson’s equations 
for singlet and triplet Green’s functions and hnally show 
how binding energy in both spin channels depends on 
inter-fermion distance r. Finally, we draw our conclu¬ 
sions in the Section IIVI and provide supplementary ma¬ 
terial in the Appendix. 


II. MODEL 

Our model system is J — J± square lattice bilayer 
Heisenberg antiferromagnet at zero temperature, where 
magnetic fluctuations are driven by interlayer coupling 
J_L (see Fig. [1]). The Hamiltonian of the undoped host 



Figure 1: Bilayer J — J± antiferromagnet model. Two black 
dotes on the top layer represent holes. 

AF reads 

= J^(s('^sW+sf^sf) + 

(iJ) 

( 1 ) 

i 

The superscripts (1), (2) in Eq. ([T]) indicate the layers, 
(i,j) denotes summation over nearest neighbour sites. 
Here = ^c\^ icr is spin of an electron at site i 

on the top plane and cJ^j^/cio -,1 is creation/annihilation 
operator of an electron with spin cr =t, j, at site f, are 
Pauli matrices. The Hamiltonian describes the antiferro¬ 
magnetic coupling in the each layer as well as between the 
two layers. It is known that without holes (half-filling) 
the model has an 0(3) magnetic QCP at J_l/J = 2.525 
(see Refs— “— ) separating the AF ordered and the mag¬ 
netically disordered phase of spin dimers. Note that since 
we consider zero temperature case, the magnetic ordering 
in the AF phase is consistent with the Mermin-Wagner 
theorem. We dope the hrst layer with two holes. For 
simplicity we set hopping integrals equal to zero, there¬ 
fore the holes are immobile. The holes interact with each 
other via magnetic fluctuations of the spins, i.e. exchang¬ 
ing by magnons. 

In the subsections A and B of the current section we 
will briefly present formalism, which describes magnon 
excitations and holc-magnon interaction on the basis of 
the bilayer model. For more detailed explanations see^^; 
a reader which is not interested in these technical details 
can go directly to the Section Hill 

A. Magnons at QCP 

Magnetic excitations in the magnetically disordered 
phase are magnons, which are also in literature called 
triplons. In the present paper we will use terms magnons 
and triplons as synonyms. To describe the magnons 
we employ the spin-bond operator mean field technique. 
This approach has been previously applied to quantum 
disordered systems such as bilayer antiferromagnets, spin 
chains, spin ladders, Hondo insulators etc— “— . It is 
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known^ii^ that this simple technique gives the position 
of the QCP at (Jj_/J)c ~ 2.31, which is close to the ex¬ 
act value {J^_/J)c = 2.525 known from Quantum Monte 
Carlo calculation a^^i^^ , series expansions^, and involved 
analytical calculations with the use of the Brueckner 
technique^^. The spin-bond technique being much sim¬ 
pler than the Brueckner technique has sufficient accuracy 
for our purposes. 

The bond-operator representation describes the system 
in a base of pairs of coupled spins on a rung, which can 
either be in a singlet or triplet (triplon) state. So, we 
define singlet and triplet operators that 

create a state at site i with spin zero and spin one, which 
is polarized along one of the axises {x,y^z). The four 
types of bosons obey the bosonic commutation relations. 
To restrict the physical states to either singlet or triplet, 
the above operators are subjected to the constraint 


by introducing infinite on-site repulsion of triplons; how¬ 
ever, this technique is quite involved. Another, more 
simple, way is to employ mean-field approach, account¬ 
ing for the constraint ([2|) via a Lagrange multiplier y in 
the Hamiltonian 

—>■ Hj j^ — -h — 1). (5) 

i 

Further analysis is straightforward. We replace sin¬ 
glet operators by numbers, (sj) = (si) = s (Bose- 
Einstein condensation of spin singlets); and diagonalize 
the quadratic in t Hamiltonian by performing the usual 
Fourier and Bogoliubov transformations 



q 


= (2) 

Oi 

In terms of these bosons, the spin operators in each layer 
and can be expressed as 

^ia ~ i ~ , (3) 

see Refi^. Substituting the bond-operator representa¬ 
tion of spins defined in Eq. ([3]) into the Hj j^ in Eq. © 
we obtain 

=Hi+H 2+H^ + Hi, 

i ^ ' 

H 2 = — {sls^^tjatja + s|sjfia4a + h.C.), 

{i,3) 

H^ = — ^ ^ -|- h.C.), 

{i,3) 

Hi = 

(h3) 

(4) 

The Hamiltonian (|T|) contains quadratic, cubic and quar- 
tic terms in magnon operators t. The most important 
for us are the quadratic terms, because they provide 
quantum criticality. The only effect due to the nonlinear 
terms H 3 and Hi is renormalization of parameters near 
the QCP, such as position of the QCP, magnon velocity, 
magnon gap and etc. This does not affect physics at the 
QCP, and therefore we will neglect these terms in further 
considerations. 

The next step for treating the Hamiltonian Q is to ac¬ 
count for the hard-core constraint ([2). It could be done 


Here N is the number of spin dimers in the lattice; the 
diagonalized Hamiltonian reads 

Hmiy^s) = 'y ) , (7) 

q 

where ujg = ^A^ and coefficients Aq = ^— 11 + 

2Js^7q, Bq = Js^"fq. Here we define 

lq = \ (cos(q„) -b cos{qy)). (8) 

The lattice spacing is set to unity. The ground state 
energy 

Eo{y,s)=Nf -i- ys'^+f3\+-'^{ujq-Aq) (9) 

k 2 q 

just shifts energy scale, and therefore is irrelevant for 
our purposes. The Bogoliubov coefficients Uq and Vq are 
given by 


Uq 



-sign(Hq) 



I 

2 ' 


( 10 ) 


The parameters y and s are determined by the sad¬ 
dle point equations: dEQ{fj,,s)/dfi = dEo{fi,s)/ds = 0. 
Solution to these equations gives position of the QCP 
at J_l/J = 2.31 and values of ’’chemical potential” 
p = —2.706 and singlet density s = 0.906. We see that 
even at the QCP s is close to unity, which again justi¬ 
fies smallness of the nonlinear terms H 3 and Hi in the 
Hamiltonian. 

The dispersion of magnons is 


Wfc = \/c2(k - Q)2 -b A2, Q = (7r,7r), (11) 

in the vicinity of the wave-vector Q, here A is the magnon 
gap and c is the velocity of magnons c = 2 Js^ = 1.64J, 
where the more precise value is c = 1.9J, see Refi^. In 
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the AF ordered phase the magnons are Goldstone bosons 
and thus necessarily gapless. On the contrary, in the 
disordered phase the gap opens up and the spin-bond 
approach gives A oc (Jj_ — J_l,c) which is not far from 
the prediction for 0(3) universality class systems A oc 
(J_L — J±,c)^ with critical index 77 = 0.71 (see Refi^). 
So, the spin-bond method provides a reasonably accurate 
description of the QCP. 

B. Hole-magnon interaction 

We dope our system with two immobile holes, by re¬ 
moving two electrons from the upper plane of the bilayer 
antiferromagnet. Hence we define the hole creation op¬ 
erator a\^ with spin projection a =t, i by its action on 
the spin singlet bond |s) 

altls) = 4t.2l0)^ altk) = CI 4 .. 2 IO) ^ ( 12 ) 

where |0) is vacuum. The electron creation/annihilation 
operator in the upper plane can be expressed in terms of 
hole creation/annihilation operators al^/uia (see Refi^), 
and after substitution in o it gives following part of the 
Hamiltonian which describes hole-magnon interaction 

Hhm = ^ ^ A (1® I ~ 

{id) 

^^z(cr 4 t] X tj]-Hcrj [t| X G]). ( 13 ) 

(®o> 

Here cTj = a\^cr The first line in the Hamilto¬ 
nian 0 corresponds to hole-magnon interaction ver¬ 
tex. The terms, describing hole-double-magnon vertex, 
which come from the second line of (fT51) will be neglected 
below, because they are irrelevant in the infrared limit. 
Performing again standard Fourier and Bogoliubov trans¬ 
formations ([S]) the Hamiltonian (1131) can be rewritten as 

Hhm ^ ( 14 ) 

i q 

The hole-magnon vertex is equal to 
J s 

9 q = --^lq{Uq+Vq). ( 15 ) 

Note, that at the QCP the vertex diverges at q —)► Q = 
(tt, tt), because of the singularity of Bogoliubov coeffi¬ 
cients Uq,Vq OC 1/—>■ OO. 

The divergence of gq is crucial for the physics of 
fermion-magnon coupling at the QCP. In fact, it re¬ 
sults in the phenomenon of spin-charge separation for 
the single-fermion problem.— The spin of the hole is dis¬ 
tributed in the power-law divergent cloud of magnons and 
in this sense is separated from the charge of the hole, lo¬ 
calized on the hole’s site. Later in the paper we will show 



Figure 2: Dyson’s equation for single hole Green’s function in 
Self-Consistent Born Approximation. Solid and waivy lines 
correspond to hole and magnon Green’s functions. 

that the infrared divergence of the fermion-magnon cou¬ 
pling at the QCP results in strong power-law attraction 
between the fermions. 

Importance of spin-charge separation for single¬ 
fermion problem at the QCP could be seen from analysis 
of analytical structure of the hole Creen’s function. Stan¬ 
dard approach in order to calculate one-fermion Creen’s 
function is to use l/Af expansion for the 0{J\f) group, 
where A/" = 3 is the number of magnon components. 
Summation of leading terms in the expansion arises in 
Self-Consistent Born Approximation (SCBA), see Fig. [21 



Figure 3: Spectral function — I/tt Im{G(e)} of a single immo¬ 
bile hole obtained in SCBA (see^). The green dashed curve 
corresponds to the QCP (A = 0), and the black solid line 
corresponds to the magnon gap A = O.IJ. Note that at the 
QCP the quasiparticle pole disappears. 

Calculations of the hole Creen’s function in the disor¬ 
dered phase in SCBA have been performed in Refs— 
The results show that away from the QCP, in the disor¬ 
dered magnetic phase, quasiparticle pole in the fermion 
Creen’s function is separated by A from incoherent part 
of the Creen’s function. But when approaching to the 
QCP the Creen’s function instead of normal pole has 
just branch cut singularity. This is a consequence of the 
infrared singularity of the coupling constant Qq. 

Spectral density of the fermion Creen’s function (see 
Fig. [3]) has inverse square root behaviour G(e) oc 
1 / V^o — £ in the vicinity of the singularity point eo and 
quasiparticle residue is approaching to zero Z cx \/A at 
the QCP. Here 

eo ~ -0.97J (16) 

is the position of the branching point of the Creen’s func- 
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tion and has a meaning of fermion energy shift due to 
interaction with magnons, where we set bare energy of 
the hole to zero. 

III. HOLE-HOLE INTERACTION, MEDIATED 
BY MAGNONS 

Now we are ready to move to the actual problem of 
magnon mediated pairing of fermions and demonstrate 
new results. Adding up magnon Hamiltonian Hm, Eq. 
© and hole-magnon interaction Hamiltonian Hhm, Eq. 
m, we arrive to effective Hamiltonian for two interact¬ 
ing holes, located at the sites with coordinates ri and 
^2, 



Figure 4: One-magnon exchange diagram that provides 
fermion-fermion interaction potential Note that 

renormalization factor y/Z should be referred to each fermion 
line. Fermion-magnon vertices also come renormalized A —>■ 
AF^^o. 


Heff 


'y ^qb\a^qa + 

<1 

E (17) 

1 = 1.2 q 




+ . . . 


The Hamiltonian (1121) is applicable only if distance be¬ 
tween holes r = |ri — r 2 | > 1, as long as we put aside 
direct exchange interaction between two neighbouring 
holes. 

The effective model m can be formulated in the lan¬ 
guage of a field theory. In fact, it is equivalent to the 
problem of two spin 1/2 fermions coupled to a vector 
field described by 0(3)-symmetric theory with La- 
grangian 

X{(f>{ri)cri + (f){r 2 )(T 2 ), 


where A is the coupling constant of fermion spin to 
magnon field. We focus only on the disordered magnetic 
phase, and therefore assume > 0. Let us note that 
in Eq. dni) we neglected self-action of magnons, hence 
term oc is dropped in Eq. (HU). 

Parameters of the Lagrangian (ITSl) could be directly 
expressed via parameters of the initial lattice Hamilto¬ 
nian ©. As an example, the coupling constant A in 
(fT51) is related to the hole-magnon vertex in the effec¬ 
tive Hamiltonian dni as Qq = Xj^2uiq. Hence, for the 
Heisenberg bilayer model A ~ 2Jy/c. This equivalence 
shows, that the problem of fermion pairing at the QCP, 
we are considering here, is generic. It has implications 
far beyond the particular bilayer model. 

In order to calculate pairing energy between two 
fermions we first consider one magnon exchange contri¬ 
bution, Fig. m According to Feynman rules we obtain 
interaction potential 


_ _ \2 ;z2p2 {cri(T2) 

KntW- ^ i-=0^2(q_Q)2^ 


A2' 


(19) 


The factor comes from for each external fermion 
line. The vertex APi^^^o comes from diagrams in Fig. [Sj 


Figure 5: The fermion-magnon vertex. 


A —>■ APtj^o- Here lu is the frequency of the exchange 
magnon, which is equal to zero. In the coordinate repre¬ 
sentation the potential reads 

Knkr) = -^^Z'^^l=oCOs{Qr){cri(T2)Ko , 

( 20 ) 

where Kq is the Macdonald function of zero’s order. The 
potential energy (r) oc ln(r) is logarithmic at small 
distances r < c/A and it exponentially decays at r > c/A 
as Vi^V(^) Spin-dependent prefactor {(T 1 CT 2 ) = 

2[S'(S' -I- I) — 3/2] is determined by the total spin of two 
fermions S and equals to —3 in singlet channel and -|-I in 
triplet channel. The potential is attractive in the state 
with total spin zero (one) when 

Pr = cos(gr) = (21) 

is negative (positive) and repulsive in the opposite case 
(r = + fy^y)- This fact has clear physical meaning 

and reflects AF character of spin correlations in the an- 
tiferromagnet. The system tends to restore AF ordering 
and the state when the spins of two interacting holes are 
aligned according to antiferromagnetic pattern (see Fig. 
HI) is energetically preferable. 

When we approach the QCP the quasiparticle residue 
as well as the magnon-hole vertex tends to zero: Z oc 
yfA —0 and T^^j^o oc A^/® —>• 0 (see discussion in Sec¬ 
tion [HB] and Refi^). Thus the single magnon exchange 
contribution given by (1^ vanishes, because the poten¬ 
tial is proportional to Z^T^^q —>• 0. Does this imply 
that pairing between fermions becomes very weak close 
to the QCP? Our answer is ”no”, on the contrary the 
pairing becomes very strong, but it is due the Casimir 
effect mechanism. 
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Figure 6: Dependence of spin channel which provides attrac¬ 
tion between holes on a mutual positioning of the holes in 
the lattice. Two holes with spins up symbolicaly represent 
triplet channel which provides negative interaction energy for 
Pr = two holes with opposite spins repre¬ 

sent singlet channel which results in attraction when Pr = —I 


An approach to evaluate Casimir interaction energy 
between two spins can be drawn from analogy with 
calculation of Van der Vaals force between two atoms. 
The most elegant way to do so was first developed by 
Dzyaloshinsky. The Van der Vaals potential is given by 
’’box diagram” with two-photon exchange between the 
atoms 1 ^ 



Figure 7: The box diagram for two-magnon exchange between 
fermions. Note that the renormalization factor \fZ should be 
referred to each external fermion line. 

We can try to follow this approach to calculate Casimir 
pairing energy between two fermions in the vicinity of the 
QCP. The diagram is shown in the Fig. 0 However, like 
in one magnon exchange diagram, we have to account for 
renormalization of external fermion lines, that results in 
additional factor 0. In such logic we again obtain 

zero for pairing energy at the QCP, and therefore we need 
another technique to solve the problem. 


A. The ’’Lamb shift” technique for calculation of 
Casimir interaction 

In this section we introduce a new technique to treat 
Casimir pairing energy. To incorporate ’’Casimir effect” 
physics, we consider composite two-fermion ’’atom”, 
which has total spin either zero (singlet state) or one 
(triplet state). Next, we calculate ’’Lamb shift” in energy 
of this composite ” atom” due to radiation of magnons as 
a function of separation between fermions. 

Let’s consider effective theory for the composite object. 


Creation operator for singlet state is 

and for triplet state 

+ «U«24,)> 


( 22 ) 


(23) 


According to the selection rules for interaction of the 
’’atom” with magnon, there are three types of tran¬ 
sitions S —>■ Ta, Ta —>■ Tp and Tq, —>• S, where 
S means singlet state and Ta denotes triplet state 
with polarization a. The only one invariant kine¬ 
matic structure that provides coupling between S and 
T, a states with emission (absorption) of one magnon is 

| 5 ST(g)^T,Q,^'s(& 9 c< + In similar way, tran¬ 
sition of the type —>■ Tp is governed by the term 

i9TT{q)£ai3-i'^T^a^T,ii{bq~f + b\^). The coefficients gsriq) 
and grriq) are coupling constants for these transitions. 
Therefore, the interaction of two-fermion system with 
magnon field in the singlet-triplet representation reads 


n = Yai3'i’T^a'^s'^gsT{q){bqp+bl^) + h.c.j + 

*£a/37^T,a^T,/3 X! 9TT{q){bq^ + ^ 57 )- (24) 

q 

The effective vertices can be calculated by evaluating 
matrix elements of the Hamiltonian (fT71) between states 

dm), dig) : 


gsriq) = grsiq) = 2 ^ 5 , sin 

grriq) = 2 c/, cos (y) ■ 

Let us define retarded Green’s function for the singlet 
and triplet state 

GT,ct/3(^2 — ^l) = —*(0|'f'T„s(i2)'I'y Q,(ti)|O)0(t2 — ^l), 

Gsit2 - h) = -^(O|d/5(<2)'I'U^l)|O)0(^2 - h), (26) 

where | 0 ) is a ground state of the system and theta- 
function is 



m 


I, t > 0; 
0, t < 0. 


(27) 


Due to the 0(3) rotational invariance the triplet Green’s 
function should be of the form GT,ap(t) = SapGTit). 
Note that our definition of the Green’s functions Gs{t 2 — 
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Ss(e) = 



ET(e) = 



Figure 8: Diagrams for singlet ^^(e) and triplet T,T{e) self¬ 
energies in SCBA. Double dashed line and dashed line rep¬ 
resent double-fermion Green’s function in singlet and triplet 
channels correspondingly. Analytical expressions for the dia¬ 
grams are given in Eq. 


ti), GT{t2 — ti) assumes that the fermions, which consti¬ 
tute the composite ” atom”, are both created at the same 
moment of time ti and then both annihilated at the mo¬ 
ment f 2 - Fourier transform of Eq. (1251) gives the Green’s 
functions in the frequency representation 

OO 

Gs,r(e) = J (28) 

0 


where ij = -|-0. Dyson’s equations for singlet and triplet 
state Green’s functions read 


Gs.T(e) = 


1 


(29) 


e - Ss,T(e) + iv' 

We use SGBA to evaluate singlet and triplet self-energy 

Ss(e) =3^! lffS'T(g)pGT(e - ujg), 
q 

^T{e) ='^\9ST{q)\'^Gsie-ujq)+ 

Q 

\9TT{q)fGT{e - ujq). 


The diagrams for the singlet and triplet self-energies are 
presented in Fig. [8l The combinator factors here come 
from contraction of the corresponding tensor structures 
of the coupling vertices in (|2^ and have a meaning of 
the number of the polarizations of intermediate state. 


B. Solution to Dyson’s equations for singlet and 
triplet states and hole-hole interaction energy 

In order to find interaction energy of two fermions, 
we numerically solve the system of two Dyson’s equa¬ 
tions (13011 in the square Brillouin zone for different inter¬ 
fermion separations r, measured in units of lattice spac- 
ings. Energy grid in our computation is Ae = 10“^ J. 
The zero approximation Green’s function is Gg^(e) = 
l/{e+iT]) and for artificial broadening we take t] = 10“^ J. 
In order to perform numerical integration in equations 
da we can directly integrate over square Brillouin zone. 


or introduce effective momentum cutoff and integrate an¬ 
alytically over angle in momentum space and then inte¬ 
grate over radial component of the momentum \q'\ = 
\q-Q\ < Aq Ri 1 numerically. We have checked that 
there is a good agreement between these two methods. 
However, the effective momentum cutoff method is much 
more efficient for numerics and provides better precision 
of the computations, therefore we mostly used the later 
approach. 

The limit r —> oo of infinite separation between the 
fermions corresponds to the case, when the vertices in 
the equations (1501) are subsituted to the averaged ones 
over q oscillations \gsT{q)?Agrriq)? -t ‘2.\gq\^ . The 
position of singularity of triplet and singlet Green’s func¬ 
tions gives us energy Eoo of the two-fermion system sep¬ 
arated by infinite distance. It is clear that in such limit 
the Green’s functions in both spin channels should coin¬ 
cide Gs(e) = Gr(e) which guarantees the same value 
for the asymptotic energy E^o in singlet and triplet 
states. So, we refer interaction energy as a difference 
Vint{r) = E{r) - Eoo. 






Figure 9: Diagrams contributing to Eoo- Top and bottom 
solid line correspond to single hole Green’s function. Diagram 
(a) is accounted in SGBA (f30l). diagram (b) is not included 
in (O and corresponds to 1/M correction to SCBA. 


We found that the value E^o ~ —1.55J (at the QPC) 
is about 20% smaller compared to doubled energy of an 
isolated single hole 2eo ~ — 1.94J, see Eq. (fTHl) . This 
difference is due to the fact that certain diagrams, which 
are presented in SCBA for single hole Green’s function, 
are not included in SCBA for the two hole Green’s func¬ 
tion Gs,t (see Fig. [5|). This deviation shows precision of 
our method, which can be improved by calculating 1/M 
corrections to SCBA dsnu- 

As it’s seen from the structure of effective vertices 

l5ST(g)P = 2^(1 - FrCOsqV), 

l5TT(g)P = 2^(1 + Pr cosq'r), 

the system’s behaviour greatly depends on the ’’parity” 
Pr = (—1)’'“:+’'!' of the inter-fermion distance. The holes 
prefer to form singlet (triplet) spin state for negative 
(positive) ’’parity” Pr at given r. 

First, let us consider the case, when the system is away 
from the QCP, A > 0. We plot spectral functions 

^S,t{c) = -Im{GsT(e)} (32) 

TT 

for singlet and triplet Green’s functions at different r, 
see FigUnKa, b). We see well defined quasiparticle peak 







3 


0.00 



Figure 10: Spectral functions ^s(e), Arie) of double-fermion 
Green’s functions in singlet and triplet channels close to the 
QCP (A = 0.08J). Panel (a) corresponds to inter-hole dis¬ 
tance r = 4 and (b) corresponds to r = 5. Blue dashed lines 
correspond to triplet state, red solid lines correspond to sin¬ 
glet state. 



in the triplet (singlet) spin channel at r = 4 (r = 5). 
However, in the opposite spin channel the peak broadens 
and submerges to continuum. This effect can be inter¬ 
preted as a formation of an excited decaying state, which 
coupled via magnons to the ground state. 

In Fig. [TT] we plot the fermion-fermion interaction 
energy versus distance. The inset displays the interac¬ 
tion energy when the system is away from the QCP (the 
magnon gap is large, A = 0.67J). Squares and triangles 
show results of our ’’Lamb-shift” technique calculations 
and solid lines represent the single magnon exchange for¬ 
mula (1^01) . There is an excellent agreement between the 
two approaches. The main part of Fig. [TT] shows the 
same quantities, but close to the QCP (the magnon gap is 
small, A = 0.08J). Here we observe a dramatic disagree¬ 
ment between the result of the ’’Lamb-shift” technique 
and the single magnon exchange potential (EUl) . The sin¬ 
gle magnon exchange approximation fails in the vicinity 
of the QCP. 

Let us now consider the most interesting case of pairing 
between fermions at the QCP (A = 0). As in the case 
of a single fermion at the QCP, the Green’s functions 
Gs{e) and Gri^) have just power-law cuts, instead of 
quasiparticle peaks, with branching point E = E{r), see 
Fig. |I2](a, b). The position E{r) of the branching point 
gives the ground state energy of the system. Spin channel 
of the ground state is specified by the spin state in which 
the Green’s function is singular at E{r). The state, in 
which the Green’s function is not singular, corresponds 
to decaying state. Imaginary part of both singlet and 
triplet Green’s functions emerges at the same branching 
point E{r) for any fixed r (see Fig. |T2)). This is due 
to transitions between the states with emission of soft 
magnons with uiq —>■ 0. In the Fig. [T2](b) we see distinct 
discontinuity of both singlet and triplet spectral functions 
at the same branching point E{r). In the Fig. |T2](a) the 
branching points are also coincide, but singlet spectral 
function has small spectral weight in the vicinity of the 
branching point. 

Our results for the interaction energy Vint{r) at the 


Figure 11: Interaction energy of two holes Vintir) at finite 
magnon gaps as a function of inter-hole distance r. Red tran- 
gles and blue squares show the results of the ’’Lamb shift” 
technique in singlet and triplet state. Red and blue solid lines 
represent theoretical prediction from one-magnon exchange 
mechanism (Unj for singlet and triplet channels. The main 
plot corresponds to small magnon gap (A = 0.08J), the inset 
corresponds to large magnon gap (A = 0.67J). 




Figure 12: Spectral functions ^^(e), Ari^e) of double-fermion 
Green’s functions in singlet and triplet channels at the QGP 
(A = 0). On the panel (a) the main plot corresponds to inter¬ 
hole distance r = 4, the inset plot to r = 20. On the panel (b) 
the main plot corresponds to r = 5, the inset plot corresponds 
to r = 21. Red solid lines show singlet state, blue dashed lines 
show triplet state. 
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QCP as a function of distance r, obtained within the 
’’Lamb shift” technique, are presented in Fig. [T31 We see 
from the data that the interaction between two fermions 
is attractive, when the ” parity” Pr is negative (positive), 
see Fig. [6] The binding becomes stronger at smaller 
inter-fermion distances r. The interaction energy has a 
power-law form 

V,nt{r) = -alr'', i/« 0.75 (33) 

with prefactor a ~ 0.3J, where v and a are found from 
the least-square fit of our numerical data. The values for 
prefactor a and power exponent v are slightly different 
for singlet (a = 0.3J, ly = 0.76) and triplet (a = 0.33J, 
ty = 0.74) cases. The variations of the values of a and ty 
are negligible within the accuracy of our calculations. 



r 


Figure 13: Interaction energy Vint{r) of two holes at the QCP, 
A = 0. Red trangles and blue squares show the results of the 
’’Lamb shift” technique for singlet and triplet state. Red and 
blue solid lines represent power-law fits Vint{r) = —ajp for 
singlet and triplet channel. The main plot corresponds to 
SCBA, in the inset we represent Vint{r), which includes first 
l/A/” correction to the SCBA. The exponent for all curves is 
approximately v « 0.75. 

In the inset in Fig. [T3]we show Vmtir) which includes 
vertex corrections to SCBA (150)) . Leading in l/AA correc¬ 
tions for the singlet and triplet self-energy S'Es and S'St 
are presented in Appendix (see diagrams in Fig. [14] and 
formulas (|A1I) , (IA2I) ). These corrections increase binding 
by about 20%, leaving critical index ly almost unchanged. 
Thus we conclude that corrections in 1/Af to SCBA do 
not change qualitative and quantitative picture, given by 
SCBA. 

From our calculations we observe very strong long- 
range attraction between fermions in the vicinity of the 
QCP. We clearly see that one magnon exchange contri¬ 
bution to the interaction energy vanishes at the QCP. 
On the contrary, accounting for multi-magnon exchange 
processes we obtain significant binding in singlet and 


triplet channels. We calculate the attraction energy due 
to multi-magnon exchange processes as a ’’Lamb shift” of 
energy of a two fermion ’’atom” due to emission of mul¬ 
tiple magnons. The fermions interact, sharing common 
’’bag” of magnetic fluctuations and reducing energy of 
fluctuations inside of the ’’bag”. Therefore, the physics 
of inter-fermion attraction in the vicinity of the QCP is 
due to ” Casimir bag” mechanism. 

IV. CONCLUSIONS 

In conclusion, we considered interaction between two 
spin 1 /2 fermions embedded in a two dimensional antifer¬ 
romagnetic system at the QCP, which separates ordered 
and disordered magnetic phases. As a model system we 
study bilayer antiferromagnet at T = 0 with two injected 
holes, in which magnetic criticality is driven by inter¬ 
layer coupling. We have shown that in the vicinity of the 
QCP the interaction between fermions can not be de¬ 
scribed by simple one-magnon exchange, unlike the case 
when the system is away from the QCP. The interac¬ 
tion mechanism is similar to Casimir effect and is due to 
multi-magnon exchange processes. To incorporate fea¬ 
tures of Casimir physics we developed a new approach, 
which we call a ’’Lamb shift” technique. We considered 
composite two-fermion ’’atom” and calculated it’s energy 
shift (’’Lamb shift”) provided by radiation of magnons. 
We found strong attraction between the fermions in spin 
singlet and triplet states depending on the ” parity” of the 
inter-fermion distance r, which is positive (negative) for 
even (odd) r. Positive (negative) ’’parity” corresponds 
to attraction in triplet (singlet) channel. The attractive 
potential has power-law form V{r) oc — 1/r^ with the 
exponent v Ri 0.75. 

We suppose that our work sheds light on the influence 
of magnetic criticality on fermion pairing mediated by 
magnons. We also believe that our results are conceptu¬ 
ally applicable to cuprates. 
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Appendix A: Leading 1/N corrections to SCBA for 
two-fermion Green’s function 
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Let us consider l/M corrections to the self-energies 
Es(e) and ST(e), calculated in Self-Consistent Born Ap¬ 
proximation (see Eqs. (l30l) L In order to do this we ac¬ 
count for vertex corrections 5Es(e) and SErie) to self¬ 
energies obtained in SCBA, corresponding diagrams are 
shown in Fig. [T31 The vertex correction to the singlet 
self-energy reads 

^S5(e) =3^ \9ST{q)\‘^\gsT{k)\'^GT,qGT,kGs,qk- 

q,k 

^'y^,9ST{q)g*sT{k)gTT{k)gTT{q)GT,qGT,kGT,qk- 

q.k 

(Al) 

The combinator factors come from contractions of the 
corresponding tensor structures of the effective vertices 
in Eq. (IMl) . In the similar way the vertex correction to 
the triplet self-energy is given by 

(5ET(e) = ^ \gsT{q)\'^\gsT{k)\'^Gs,qGs,kGT,qk- 

q.k 

2 ^ gTTiq)gTTik)9sTiq)9ST{k)GT,qGT,kGs,qk+ ^^2) 

q,k 

2 ^ \9TT{q)\^\gTT{k)\'^GT,qGT,kGT,qk, - 

q.k 
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